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ABSTRACT 



An optimal integral control design program was applied 
to a linearized state variable model of the FFG-7 ship class 
gas turbine and Controllable Reversible Pitch (CR?) propeller 
main propulsion system. Various combinations of output 
parameters were investigated in an attempt to produce a 
feasible control design. Only one acceptable design was 
achieved which did not violate any physical constraints. 
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I. 



INTRODUCTION 



A. OBJECTIVE 

Advanced technology propulsion systems require a control 
system design that is able to achieve the desired performance 
at the greatest possible efficiency. One particular scheme 
that shows great promise in this regard is the optimum integral 
control design program, CONSYN, developed in Ref. 1. 

This study was undertaken to examine the validity and 
effectiveness of applying CONSYN to the FFG-7 ship class gas 
turbine propulsion system model obtained from Ref. 2. 



B. COMPUTER FACILITIES AND SOFTWARE UTILIZED 

The above objective was implemented by simulation on the 
IBM 360/67 computer of the W. R. Church Computer Center at 
the Naval Postgraduate School. The following programming 
software was required: 

Various routines in the International Mathematical and 
Statistical Library (IMSL) 

The IBM developed Continuous System Modeling Program 
(CSMP) 

’ The Control Program for Engineering Synthesis with 
Constrained Function Minimization (COPES/CONMIN) de- 
veloped by Garret N. Vanderplaats of the NASA Ames 
Research Center 
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The routines developed for CONSYN in Ref. 1 as modified 
by Ref. 3 and this study, are listed in Appendix A. 
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II. PROPULSION PLANT MODEL 



The propulsion system for the OLIVER HAZARD PERRY (FFG-7) 
class patrol frigates consists of two General Electric 
Corporation LM2500 Marine Gas Turbines, coupled through a 
reduction gear to a controllable reversible pitch (CRP) 
propeller manufactured by the 3ird-Johnson Company. Since 
the gas turbine is a highly non-linear device, the linearized 
state variable model developed in Ref. 2 is valid only for 
small changes from a specified steady state operating point. 
This model was obtained from a coupled ship and propulsion 
plant dynamic non-linear model consisting of the following 
features : 

A non-linear model of the LM2500 engine developed by 
General Electric. 

A non-linear model of the ship, propeller and propul- 
sion train provided by the Bath Iron Works Company . 
Secondary flow losses considered. 

Bleed flow for extraction and turbine cooling considered. 
Inlet and exhaust losses incorporated. 

Pressure and temperature dynamics neglected. 

* Nested loop balancing utilized to obtain various steady 
state operating points. 

In matrix equation form, this linearized state variable 
model of the various engine and ship parameters is expressed 
as : 
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X 



AX + BU 



Y = CX + DU 
where 

X is the system state vector, and consists of three 
components : 



1 — 1 
X 


= NG - NG 

ss 


; NG — Ga.S 


Generator Speed 


(rpm) 


X 2 


= NPT - NPT 

ss 


; NPT = Power Turbine Speed 


(rpm) 


x 3 


= V - V 

ss 


; V = Ship 


Speed (knots) 





U is the system input vector, and consists of four 
components : 



u, = WF - WF 
1 ss 



WF = Fuel Flow (lbm/hr) 

3 = Stator Guide Vane Angle (deg) 
WB = Compressor Bleed Flow (lbm/s) 
P/D = Propeller Pitch/Diam Ratio ( 
Y is the system output vector, and consists of six 
components : 



u~ = 3 - 3 cc: 

2 . s s 

U-, = WB - WB 
3 SS 

u 4 = P/D - P/D ss 



y, = SMC = SMC 



SS 



y- = T4 - T4 



ss 



y 3 = QPT - QPT ss 

y. = P3 - P3 
1 4 ss 

y 5 = HPPT - HPPT 
= T5 - T5 



ss 



SMC = Compressor Surge Margin (%) 

T4 = Combustor Outlet Temp (°R) 

QPT = Power Turbine Torque (ft-lbf) 

P3 = Compressor Total Disch Press (psia) 
HPPT = Power Turbine Power (hp) 

T5 = Power Turbine Inlet Temp (°R) 



1 6 ss 

and the ss subscripts indicate values at a specific steady state 
operating point, the same point used for determination of the 
constant A, B, C, D matrix elements. 

This study considered single engine up-power transients of 
70% to 75%, and 80% to 85% of maximum gas generator speed. In 
both cases, the steady state operating point selected was that 
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for the higher power, or final condition. The necessary 
state variable data to accomplish this study was obtained 
from computer printouts generated in conjunction with Ref. 2, 
and is contained in Appendix B. Physical constraints on 
inputs and input rates are also contained in Appendix 3. 
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III. CONTROL SYSTEM DESIGN WITH CONSYN 



A. BASIC PHILOSOPHY 

The CONSYN (Control Synthesis) design program presented 
in Ref. 1 is a further refinement of linear regulator theory 
developed in Ref. 4. It is based on the initial assumption 
that the plant dynamics can be represented by the linearized 
state variable matrix equations 

X = AX + BU 

Y = CX + DU 

where, as previously discussed, all states are perturbed 
from their steady state operating points. This is illus- 
trated in block diagram form in Fig. la. This original plant 
is then augmented with integrators to achieve open loop inte- 
gral control as shown in Fig. lb. Application of linear 
quadratic regulator theory to this augmented system leads 
to the minimization of a performance or cost index, J, 
defined by the matrix equation 

fc f 

J = f (Y T QY + U T RU) dt 
o 

where Q is a constant, diagonal, positive semi -definite 
matrix used to weight output deviation from final values, and 
R is a constant, diagonal, positive definite matrix used to 
weight input rate deviations from final values. 
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The linear regulator problem solution yields an optimal 
state feedback regulator which follows the matrix equation 
control law 

U = Z - G 1 X - G 2 U 

and is illustrated in Fig. lc , where Z is the demand vector. 
Finally, in order to achieve integral control commands gen- 
erated by the error between demand and output, the system 
shown in Fig. lc is converted to that shown in Fig. Id. This 
final system configuration has the same time response as that 
shown in Fig. lc but has the added advantage of assuring zero 
steady state error even if system degradation alters the 
elements of coefficient matrices. 

Since that choice of Q and R matrix elements for the 
solution of the linear regulator problem is arbitrary, 
COPES/CONMIN logic is then utilized to find the optimum Q and 
R which minimizes the cost index J without violating any of the 
imposed physical constraints. 

In summary then, the CONSYN design program applies a 
linear regulator design form to the given output control 
problem using the augmented state vector 

X* = [-■§-] . 

Block diagram methods along with suitably partitioned matrices 
are utilized to arrive at the desired structure of the optimal 
control system. 
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B. PROGRAM FLOW LOGIC 



The CONSYN design program operates in three distinct 
segments: Initial Design, Optimization, and Final Design 

Output. A summary of the significant logic and operations 
within each of these segments follows. 

1 . Initial Design 

a. Input the following: 

(1) A, B, C, D matrices 

(2) Initial guess for appropriate Q and R matrices 

(3) Initial condition X and U vectors, which are 
then combined to form the augmented state 
vector X* = [-—] 

(4) System constraints 

b. Calculate system open loop characteristics 

c. Calculate P, the solution to the steady state 

matrix Riccati equation" 

A* P + PA* - PB* R~ 1 3* P + C* Q C* = 0 

where A* = [-5 j§-] B* = [-2-] c* = [C j Dj 

d. Calculate G the optimal state feedback gain matrix 

G = R -1 3* P = [G 1 |g 2 ] 

■^The detailed development of this Riccati equation and 
the following G, H, and L gain relationships in the partitioned 
matrix form is found in Refs. 1 and 4. 
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e. Calculate H and L the integral control gain matrices 

[l|h] = G [-Ai-A] _1 

I 

where E = [i !o] and is required so that the number 

i 

of outputs fed back for comparison with the demand 
vector is equal to the number of inputs. 

f. Simulate time response of system. 

g. Print results of all of the above. 

2 . Optimization 

Vary Q and R matrix elements under the direction of 
COPES/CONMIN logic to achieve the smallest cost function 
value without violating any constraints. 

3 . Final Design Output 

Print complete results of system charac teris tics 
associated with the optimum Q and R matrix elements. 
Since the CONSYN program does not contain a plotting routine 
for system time response simulation, a separate program must be 
utilized. For this study, the IBM CSMP in conjunction with a 
VERSATEC plotter was employed. A derivation of the equations 
required for this simulation is contained in Appendix C. 



C. MODIFICATIONS TO THE ORIGINAL PROGRAM 

As discussed in Ref. 3, the following two significant modi- 
fications were made to the original CONSYN program: 

The Kleinman technique for the solution of the matrix 
Riccati equation in subroutine OPGAIN was replaced by an 
eigenvalue method. 
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* The fixed step integration solution in subroutine PEAK 
for system simulation was replaced by a variable step 
integration routine. 

In using CONSYN with these changes for this study, diffi- 
culty was encountered with the variable step integration routine. 
Numerical underflows, overflows, and step size decimation were 
encountered to varying degrees in all attempts to use it. Con- 
sequently, the fixed step integration method was returned to 
subroutine PEAK, along with a logic modification to determine 
and print out the number of terms required for the matrix 
series convergence of PHI = e . Here, the matrix F is 
defined as 

F = A* - B* G 

where the matrices A*, B* and G are as previously defined. 
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Fig. la. Original Plant 
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Fig. lb. Open Loop Integral Control 



Figure 1. Integral Control Synthesis 



18 



± + 



i 




y 



Fig. lc . Closed Loop Optimal Integral Control 




Fig. Id. Closed Loop Optimal Output Integral Control 



Figure 1. (Cont'd) 



19 



IV. RESULTS 



A. GENERAL 

The overall results were disappointing. Of all the 
basic design forms and variations evaluated, only one was 
found which did not violate one or more physical constraints . 
In addition, the study was hampered by the relatively long 
computer execution times required for the optimization of 
a single proposed design (approximately 20 minutes) , and the 
excessively long computer turn-around times (approximately 
10 hours) . In all, at least 150 difference designs were 
attempted, with a total computer execution time expenditure 
of over 50 hours. A summary of the design forms and varia- 
tions considered is given in Tables I and II. 

Initial efforts to find an acceptable design were con- 
centrated on design form numbers 2 and 4 in Table I. These 
forms were selected because they appeared to be the most 
analogous to the design evaluated in Ref. 5 for the F401 
turbofan jet aircraft engine with respect to those outputs 
designated for feedback.^" All proposed designs of these 
forms consisted of transients between 70% and 75% of maximum 
gas generator speed, variable bleed, and no state variable 
scaling. All variations considered resulted in violation 



■^See Appendix D for a more detailed discussion. 
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of one or more physical constraints. Examples most often 
noted were: 

Maximum propeller pitch exceeded 

Maxim-urn propeller pitch rate exceeded, either in the 
positive or negative direction 

Negative bleed 

Design form number 4 was then exercised at various 
transients between 30% and 85% of maximum gas generator 
speed with similar types of physical constraint violations. 

In an effort to determine why negative bleed was some- 
times occurring, G, J. Michael, one of the authors of Refs. 

4 and 5 was contacted. The following information and guidance 
was obtained from this communication: 

Although variable bleed was used as a control in Ref. 4, 
it also sometimes assumed negative values due to mathe- 
matical anomalies in the engine model that could not be 
eliminated. As a result, bleed was not used as a 
control input in Ref. 5. 

Although the primary motivation for state variable 
scaling in Refs. 4 and 5 was to avoid publishing 
proprietary information in the open literature, doing 
so in the case of this study might eliminate additional 
computational anomalities. 
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TABLE I 



DESIGN 


FORMS 


EVALUATED 










Basic Form 


I .D. 


Number 




2 


4 


5 


6 


7 


3 


U 1 


U 1 


U 1 


U 1 


U 1 


U 1 


U = U 2 


U 2 


u 2 


u 2 


u 2 


u 2 


U 3 


U 3 


U 4 


U 4 


U 4 


U 4 


U 4 


U 4 










u 2 


u 2 


U 2 


U 2 


u 2 


U 1 


u 4 


X 3 


X 3 


y 3 


V - 
^ 0 


u 4 


y 3 


y 3 


y 6 


y 6 


y 6 


x 3 


i 

1 

U>l 
>1 1 

1 

II 


y 6 


y 3 


X 3 


X 3 


y l 


y i 


y l 


y l 


y l 


y l 


y 2 


y 2 


y 2 


y 2 


y 2 


y 2 


v. 


y 4 


y 4 


y 4 


y 4 


y 4 


y 5 


y 5 


y 5 


y 5 


y 5 


y 3 


y 6 



NOTE: The first n elements of the Y vector, where n is the 

number of inputs, are those that are to be fed back 
for comparison with the demand vector . 
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TABLE II 



DESIGN FORM VARIATIONS EVALUATED 



1. Initial to final Power Level, % of Maximum Gas Generator 
Speed : 

a. 70 - 75% 

b. 80 - 85% 

2 . Bleed : 

a. Variable 

b. Constant 

3. Settling time - 20 to 200 seconds. 

4. Ship speed constraint with respect to settling time: 

a. Within + 2% of steady state value 

b. Not constrained 

5. State variable scaling or normalization with respect to 
the difference between initial and final steady state 
values : 



a. Scaled 

b. Unsealed 

6. Elements of the Q and R performance index weighting 
matrices allowed to be changed by the optimization 
routine : 

a. All Q and R elements 

b. r^ (the weight of fuel flow rate error) fixed at 
1.0, all other Q and R elements allowed to change. 
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Design form numbers 5, 6, 7 and 8 were then exercised 
at various transients between 80% and 85% of maximum gas 
generator speed, constant bleed, and with state variable 
scaling.'*' The following results were noted: 

An acceptable design which did not violate any of the 
input or input rate constraints was achieved with 
design form number 6, with the performance index 
weight on fuel flow rate error fixed at 1.0. 

There was no significant difference in system response 
obtained from design form numbers 5, 6 and 7. 

For those designs which did not impose a specif ication 
on ship speed with respect to requested settling time, 
the optimization routine reduced the degree of initially 
violated constraints during the first several intera- 
tions . The remaining iterations produced no further 
changes in the degree of constraint violation or in 
the value of the objective function. 

For those designs which did impose a specification on 
ship speed with respect to requested settling time, 
the degree of initially violated constraints was 
reduced at the expense of violating another constraint 
that was not initially violated. 



See Appendix E for the derivation of the scaled state 
variable data with constant bleed. 
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The number of terms required for the series convergence 
Ft • 

of PHI = e in subroutine PEAK was much less than 
before state variable scaling was utilized; usually 
fewer than ten terms were needed for convergence. 

Deviations from actual values were noted for the 
linearized model calculations for outputs at the initial 
condition. The worst occurred for power turbine power 
which assumed a negative value. 



B. SPECIFIC 

A summary of the design data associated with the one 
acceptable design achieved is presented in Table III. Plots 
of system time response simulation for this design are shown 
in Figures 2a through 2f. Plots of state open loop response 
to a step input corresponding to the final steady state input 
values are shown in Figures 3a and 3b. It is noted that with 
the optimum design, all engine parameters reached their final 
steady state values within 6 seconds, and the ship speed 
achieved 99% of its final steady state value after 100 seconds. 

Summaries of performance data for various selected designs 
are presented in Table IV, with the degree of constraint 
violation indicated. 

A comparison of actual output values with those predicted 
by the linearized model at the initial condition is given in 
Table V. 
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TABLE III 



OPTIMAL CONTROLLER DATA, DESIGN 6-2F 



Performance Index Weighting Factors 



Output Error 


Input 


Rate Error 


*1 


q 2 q 3 


q 4 q 5 q 6 


q 7 


q 8 


r l 


r 2 r 3 


(u 2 ) 


(y 3 ) (y 6 ) 


(x 3 ) (y (y 9 ) 


(y 4 ) 


(y 5 ) 


(u ] _) 


{u 9 ) (u 4 ) 


.66 


.37 .26 


7x10 1 6x10 ^ .56 


.72 


. 023 


1.0 


.32 .28 


The 


H Matrix 














.3032 


.5066 .9441 












1.533 


.03452 -.02609 












-.1797 


.9716 -.8782 











The L Matrix 

-.3585 -.02830 

.1184 -.03476 

.8133 -.3113 



-.0006817 
- .001542 
-.2850 



NOTE: The complete computer output of the optimization run for 

this design is contained in Appendix F. 
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TABLE IV 



PERFORMANCE DATA SUMMARY FOR VARIOUS DESIGNS 
Design Code 



Constraint 


Limit 


2-lA 


< 
t — \ 

I 


4-1B 


4-2B 


WF 

max 


7600 


1640 


1640 


1350 


33 60 


WF . 
min 


900 


1320 


1320 


1320 


2360 


s 

max 


40 


33.3 


33.3 


33.3 


23.4 


^min 


0 


29.0 


27 .5 


29.4 


13.0 


WB 

max 




• — i 
'X) 

• — i 


2.54 


2.03 


10.4 



P/D 
' max 


1.43 


(T05) 


€ 3 )( 0 °) 


1, 


.43 


P/D . 
min 


.01 


1.29 


1.29 


.93 


1 . 


.17 


WF 


+ 3200 


30 . 5 


34.7 


1090 


1 ? 


520 


3 


+ 59 


-45.3 


-43.3 


-10.2 


- i ; 


1.1 


WB 




11.0 


19.1 


2.53 


24 6 


P/D 


+ .12 


(+ 59 ) 


(fT 02) 


(2730) ( 


_2 , 


_9f; 



6-2C 


6- 2D 


6-2E 


6-2F 


3-2E 


3 -2F 


3370 


4170 


3360 


3360 


3950 


3910 


23 60 


2360 


2360 


2360 


2360 


2 3 6 0 


23 .4 


23 . 4 


23.4 


23.4 


23.4 


23.4 


13.3 


15.0 


13.2 


13.2 


11.3 


13.1 








n 








t ul I O ^ u 1 1 L — • 












n 








C O li S u Q li L * 








1.43 


1.4 ’ 


1.43 


1.43 


1.43 


1.43 


1.05 


.30 


1.15 


1.17 


1.15 


1.12 


611 


(4t7o) 


(4650) 


1330 


(3sT q) 


3130 


-12.7 


-1.3 


-17 .2 


-16 . 5 


-3.7 


-6.2 



constant: = 0 

(£Zza> -- 12 G 5 CD) 




NOTE: Circled items indicate constraint limit exceeded. The design code 

used above is explained as follows : 

(1st Digit) - (2nd Digit) (letter) 

i t | 

Design form Power Code Variation Code 

from Table I 1 = 70-75% from next page 

2 = 30-35% 
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TABLE IV (Cont'd) 



Variation Code (Also see Table II) 



A B C D E F 



Bleed 



Variable 

Constant 



X X 



X X X X 



State Variables 
Scaled 
Unsealed 



X X X 



XXX 



Q and R Weights 
Variable 
r^ Fixed 

Ship Speed 

Constrained 

Unconstrained 



X X X X X 



X 



XXX 

X X X 
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TABLE V 



COMPARISON OF INITIAL CONDITION ACTUAL OUTPUTS WITH THOSE 
COMPUTED BY THE LINEAR MODEL 



Output 


Actual Value"'' 


Computed Value 


SMC 


32.39 


30 .00 


T4 


1776 


1853 


QPT 


9706 


8132 


P3 


98 . 53 


91 . 09 


HPPT 


3219 


-1795 


T5 


1316 


1365 


'^'Obtained 


from data generated in 


conjunction with Ref. 2. 
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Figure 2a. Gas Generator Speed and Power Turbine Speed 
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Figure 2b, Ship Speed and Propeller Pitch to Diameter Ratio 
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Figure 2c. Fuel Flow Rate and Stator Guide Vane Angle 
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Figure 2d. Combustor Outlet Temperature and Power Turbine Inlet Temperature 
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Figure 2e. Compressor Surge Margin and Total Discharge Pressure 
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Figure 2f. Power Turbine Torque and 
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Figure 3a, Gas Genera Lor Speed and Power Turbine Speed 
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Figure 3b. Ship Speed 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. Some form of state variable scaling is necessary to 
improve the" efficiency of the numerical calculations. 
Scaling causes the A, B, C, D matrix elements to be 
much closer to the same order of magnitude. 

2. If bleed is to be utilized as a control variable, 
the failure of the mathematical model to prevent it 
from assuming negative values must be corrected. 

3. Although the system's non-linearity was in evidence 
at the initial condition state as shown in Table V, 
the mathematical model became a much more accurate 
representation of the actual system as the final 
steady state operating point was approached. 

4 . Optimal Controller design with the CONSYN program is 
by no means a completely automated (or cookbook) 
process. Much manual interaction with options and 
variable parameters must be employed in order to 
achieve a feasible design. 

5. Input and input rate constraints were design limit- 
ing. As long as these constraints were not violated, 
there were no unreasonable excursions of any state or 
output variables. 
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B . RECOMMENDATIONS 



1. The specific reason (s) why the optimization routine 

is so frequently unable to reduce violated constraints 
beyond a certain point must be investigated and 
corrected before any practical controller design by 
this method can be achieved. 
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APPENDIX A 





LIST OF CONSYN SUBPROGRAMS 


Name 


Purpose 


ANAL I Z 


Integral control design and analysis (pre 
viously named CYCLE in Ref. 1) 


CONTL 


Controllability Check 


COST 


Evaluation of Performance index and deter 
mination of maximum eigenvalue 


DET 


Matrix determinant 


INVERT 


Matrix inverse 


LIAP 


Solution to bilinear symmetric Lyapunov 
equations 


MULT 


Matrix multiplication 


MULTV 


Vector multiplication 


MULTZ 


Complex matrix multiplication 


OBSERV 


Observability check 


OPGAIN 


Solution to optimal linear regulator 


PEAK 


System simulation 


PHIOFT 


State transition matrix 


POLY 


Coefficients of characteristic polynomial 


PRINTM 


Matrix print, double precision 


PRINTN 


Matrix print, single precision 


RANK 


Matrix rank 


ROOT 


Polynomial roots 
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APPENDIX B 



PROPULSION PLANT DATA 

1. State Variable Representation at 70%, 75%, 
80%, and 85% of Maximum Gas Generator Speed 





70% 


75% 


80% 


85% 


States 


NG (rpm) 


6879 . 


7370. 


7862. 


8353 . 


NPT (rpm) 


1058 . 


1313 . 


1742 . 


2300 . 


V (knots) 


9.688 


11.91 


15.64 


20.40 


Inputs 


WF (lbm/hr) 


1319 . 


1640 . 


2361 . 


3857 . 


8 (deg) 


33.29 


30.53 


23.37 


13 . 20 


WB (lbm/S ) 


.02 


. 02 


.02 


.02 


P/D ( * * * ) 


1.43 


1.43 


1.43 


1.43 


Outputs 


SMC (%) 


32.63 


33.21 


32.39 


27.76 


T4 (°R) 


1613 . 


1663 . 


1776 . 


2002 . 


QPT ( ft-lbf ) 


3745. 


5627 . 


9706 . 


17040 


P3 (PSIA) 


54 . 02 


69.24 


98 . 53 


144 . 6 


HPPT (hp) 


754.4 


1407 . 


3219. 


7461. 


T5 (°R) 


1284 . 


1277 . 


1316 . 


1456 . 
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APPENDIX B (Cont'd) 



2. A, B, C, D MATRIX ELEMENTS FOR 75% 



OF 


MAXIMUM GAS 


GENERATOR 


SPEED 


The A Matrix 


-.2709 


-.1186 


0 . 0 




.1295 


-.9364 


57.71 




0 . 0 


.0002429 


-.02793 




The B Matrix 


.7287 


37.74 


-28 . 29 


0.0 


.08656 


-8.978 


-7.861 


-526. 


0.0 


0.0 


0.0 


.1114 


The C Matrix 


. 01078 


-.0003271 


0 . 0 




-.2160 


.0008184 


0 . 0 




3.171 


-2.473 


0 . 0 




. 02112 


.0001619 


0 . 0 




.7928 


.4526 


0 . 0 




-.2289 


.01866 


0 . 0 




The D Matrix 


-.01811 


-.5678 


2.453 


0.0 


.4905 


9.434 


10.21 


0.0 


2.120 


-219.9 


-192.5 


0 . 0 


.008820 


-1.405 


-1.195 


0 . 0 


.5300 


-54 .97 


-48.13 


0.0 


. 3635 


11.49 


11.85 


0.0 
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APPENDIX B (Cont'd) 



3. A, B, C, D MATRIX ELEMENTS FOR 35% 
OF MAXIMUM GAS GENERATOR SPEED 



The A Matrix 


-1 . 212 


-.1923 


.3466 


-1.711 


o 

• 

o 


.000429 


The B Matrix 


.6284 


68.15 


.1194 


-14.23 


O 

• 

o 


o 

o 


The C Matrix 


.01097 


0 . 0 


-.2715 


0 . 0 


8.489 


-5.786 


.04678 


0.0 


3.717 


7.101 


-.2484 


.01736 


The D Matrix 


-.007521 


-.3569 


.2423 


7.875 


2.925 


-348.4 


.008168 


-1.802 


1.281 


-152.6 


.1855 


7.277 



0.0 
100 . 3 
0 -.05131 



-26.72 


0.0 


-10 .79 


-1634 


0.0 


. 344 5 



0.0 




0 . 0 




0.0 




0 . 0 




0 . 0 




o 

• 

o 




1.229 


0 . 0 


7.273 


0.0 


-264.3 


o 

• 

o 


-1.349 


0 . 0 


-115.8 


o 

• 

o 


6.693 


0 . 0 
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APPENDIX B (Cont’d) 



4 . SCALED A, B, C, D MATRIX ELEMENTS FOR 35% 
MAXIMUM GAS GENERATOR SPEED, CONSTANT BLEED 



The Scaled 


A 


Ma tr ix 




- 1.212 




-.2185 


0 . 0 


.305 




-1.711 


.3556 


O 

• 

o 




.05029 


-.05181 


The Scaled 


B 


Matrix 




1.915 




1.412 


0.0 


.3201 




-.2594 


-2 .928 


0 . 0 




o 

• 

o 


.07237 


The Scaled 


C 


Matrix 




1.163 




O 

• 

o 


O 

• 

o 


-.5899 




0 . 0 


0 . 0 


.5683 




-.4402 


0 . 0 


.4986 




0 . 0 


O 

♦ 

O 


.4302 




.9341 


0 . 0 


-.8712 




. 06919 


0 . 0 


The Scaled 


D 


Matrix 




-2.43 




-.7839 


0 . 0 


1.604 




.3544 


o 

• 

o 


.5966 




-.4831 


0.0 


.2652 




-.3978 


o 

• 

o 


.4518 




-.3659 


O 

• 

o 


1.982 




.5286 


o 

• 

o 
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APPENDIX B (Cont'd) 



5. PHYSICAL CONSTRAINTS 

Inputs 

900 < WF < 7600 lbm/hr 
0 < 3 < 40° 

.01 < P/D < 1.43 

Input Rates 

-3200 < WF < 3200 lbm/hr/S 
-59 < 3 < 59 °/S 
-.12 < P/D < .12 1/S 



Note 



( 1 ) 

( 1 ) 

( 2 ) 



(3) 

( 3 ) 

(4) 



NOTES: (1) 

( 2 ) 

(3) 

(4) 



Approx, range between idle and full power 
condition . 

Model does not allow for negative pitch. 
Obtained from G.S. 

Computed from data provided by Bird Johnson 
Company as follows : 



Propeller Diameter = 16.5 ft 

Maximum Ahead Pitch = 23 . 5 ft 

Maximum Reverse Pitch = 14.7 ft 

Minimum Time for Pitch Change from 
Full Reverse to Full Ahead = 20 sec 






23 .5 ,_14 .7, 

16.5 1 16.5 
20 S 



.12 1/S 
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APPENDIX C 



DERIVATION OF EQUATIONS USED 
FOR SYSTEM SIMULATION WITH CSMP 

The following matrix equations may be obtained from 
inspection of the system block diagram in Fig. id: 

X = AX + BU 
Y = CX + DU 
U = /[ H (Z-EY) ]dt - LX 

or, upon differentiating, 

U = H(Z-EY) - LX 

However, the demand vector Z is equal to zero since all 
desired output deviations are zero at the final steady 
state operating point which has been defined as the 
origin . Thus , 

U = -(HEY + LX) 

Recalling that the matrix E was defined as 
E = [I |0], 

the above matrix equations can, therefore, be expanded 
into the following form for the 3-state, 3-input, and 
8-output system considered: 



X 1 


a ll 


X 1 




a 12 


X 2 


4- 


a 13 


x 3 


4- 


b n 


U 1 


4* 


b 12 


u 2 


+ °13 


U-, 


X 3 


= a 31 


X 1 


4* 


a 32 


X 2 


4* 


a 33 


x 3 


+ 


b 31 


U 1 


+b 3 2 ' 


U 2 


- b 33 


u ^ 

Jj 


y l 


= “11 


X 1 


+ 


C 12 


x 2 


4* 


c 13 


X 3 




d ll 


U 1 


+ 


d 12 


u 2 


+ d 13 


U 3 


00 


C 81 


X 1 


4* 


C 82 


x 2 


4* 


C 8 3 


x 3 


4- 


d 81 


U 1 


4- 


d 8 2 


U 2 


+ a 83 


U 3 



46 



APPENDIX C (Cont'd) 



(hu y L + h 12 

(h 31 Y 1 + h 32 



+ h 13 


y 3 


+ 1 L1 x 1 


+ L I2 


X 2 


+ X 13 


x 3 ) 


+ h 33 


y 3 


+ 1 31 x L 


+ 1 3 2 


x 2 


+ 1 3 3 


i 3 ) 
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APPENDIX D 



RATIONALE FOR DESIGN FORM SELECTION 

In Ref. 5, the following outputs were selected for 
feedback in the case of the F401 turbofan jet engine 
controller design: 

1. Fan inlet guide vane angle (FIGV) 

2. Rear compressor variable vane angle (ROW) 

3. Thrust 

4 . High turbine inlet temperature (HTIT) 

Those outputs selected for feedback in design form number 
2 of this study were: 

1. Propeller pitch to diameter ratio (P/D) 

2. Compressor stator guide vane angle ( 3 ) 

3. Power turbine torque (QPT) 

4 . Power turbine inlet temperature (T5) 

A comparison of these designs follows: 

1. The marine propulsion gas turbine engine, of 
course, has no turbofan. However, the ?/D was 
considered analogous to the FIGV since both 
directly relate to the loading of the output 
turbine. In the aircraft case, the output 
turbine drives the fan; in the marine case, 
the output turbine drives the propeller. 

2. 3 is completely analogous to RCW since this is 
the compressor vane control in both cases. 
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APPENDIX D (Cont’d) 



3. The thrust in the aircraft case is a measure of 
the total power produced. In the marine case, 
either ship speed, power turbine power, or power 
turbine torque is analogous. Torque was selected 
since it was considered the easiest of the three 
to measure accurately. 

4. HTIT is ' equivalent to combustor outlet temperature. 
However, the high temperature involved shortens 
transducer lifetime, and is variable around the 
perimeter of the combustor. To eliminate these 
difficulties, power turbine inlet temperature, 
instead, was chosen for feedback. 
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APPENDIX E 



DERIVATION OF SCALED STATE VARIABLE DATA 
WITH CONSTANT BLEED 



From the matrix equation, 

X = AX + BU 

consider one of the scalar equations 

X I = a ll X 1 + a 12 X 2 + a 13 X 3 + b ll U 1 + b 12 U 2 

+ o 13 u 3 + b 14 u 4 

f s 1 s 

Since the x and u are deviations from the steady 
state operating point, and bleed is constant, then u 3 = 0. 

Now define the following scaled variables: 
x i u i 

i Ax. l Au. 4 4 

r l 

* s 

where the A refer to the difference between initial and 
final conditions, and u^ will not be scaled since its 
initial and final values are equal. 

Then 

Ax. x. = a.. Ax. x. + a. 0 Ax_ x 0 + a. Ax n x, 

1 1 11 1 1 12 2 2 13 3 3 

+ b lx Au x u x + b 12 Au 2 (i 2 + b 14 u 4 

Dividing through by Ax^ results in 

• Ax 2 ^ Ax 3 ^ Au-, ^ 

*1 - a n *1 + a 12 ixT *2 + a i3 53V x 3 + b ll 5337 U 1 



iu 2 . b 14 . 

12 Ax^ U 2 Ax^ J 4 
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APPENDIX E (Cont'd) 



Carrying this procedure out with the other state and 
output equations results in the scaled matrix equations 



X 

Y 

where 



A 



AX + BU 

A A A A 

CX + DU 



a 



11 



a 21 






Ax 9 


Ax. 


a, - — - 
12 Ax x 


a, -3 ^ — 
13 Ax 




Ax 


a 22 


a 23 Ax 



a 3 2 Ax-, a 33 





b ll 


Au^ 

Ax-^ 


b 12 


Au- 

Ax l 


-r 

i — 1 

X) 


1 

Axi 


3 = 


b 21 


Au^ 

Ax 2 


b 22 


Au-, 

Ax 2 


b 24 


1 

Ax 2 




b 31 


Au^ 


b 32 


Au 2 




1 




Ax^ 


Ax 2 


b 34 


Ax 3 



51 



APPENDIX E (Cont'd) 





i — 1 

<3 




AX 2 




Ax^ 


C 11 


AYl 


C 12 


ayi 


C 13 


A Yl 




Ax-^ 




Ax 2 




Ax^ 


C 21 


Ay 2 


C 22 


Ay 2 


C 23 


Ay 2 




Ax, 




Ax 0 




Ax_ 


C 61 


**6 


C 62 


A ^6 


C 63 






Au^ 




Au 2 




1 


i — 1 
i — 1 

3 


Ay 1 


d 12 


A Yl 


d 14 


A Yl 


d 21 


Au-^ 

Ay 2 


d 22 


A u 2 

Ay 2 


d 24 


1 

a y 2 


d 


Au-^ 




Au 2 


d 


1 


61 


A ^6 


62 


A ^6 


64 


A ^6 
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APPENDIX F 



COMPUTER OUTPUT OF OPTIMIZATION RUN FOR DESIGN 6-2F 
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DESIGN VARIABLES 
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SIMULATION RUN TIME *= /2 .000 
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PEAK VALUES T 1 ME OF OCCURR ANCE 
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